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Abstract 

Spinor ultracold gases in one dimension represent an interesting example of strongly correlated 
quantum fluids. They have a rich phase diagram and exhibit a variety of quantum phase transitions. 
We consider a one-dimensional spinor gas of bosons with a large spin S. A particular example is the 
gas of chromium atoms (5 = 3), where the dipolar collisions efficiently change the magnetization 
and make the system sensitive to the linear Zeeman effect. We argue that in one dimension the 
most interesting effects come from the pairing interaction. If this interaction is negative, it gives 
rise to a (quasi) condensate of singlet bosonic pairs with an algebraic order at zero temperature, 
and for {2S + 1) ^ 1 the saddle point approximation leads to physically transparent results. Since 
in one dimension one needs a finite energy to destroy a pair, the spectrum of spin excitations has 
a gap. Hence, in the absence of magnetic field there is only one gapless mode corresponding to 
phase fluctuations of the pair quasicondensate. Once the magnetic field exceeds the gap another 
condensate emerges, namely the quasicondensate of unpaired bosons with spins aligned along the 
magnetic field. The spectrum then contains two gapless modes corresponding to the singlet-paired 
and spin-aligned unpaired bose-condensed particles, respectively. At T=0 the corresponding phase 
transition is of the commensurate-incommensurate type. 

PACS numbers: 05.30. Jp, 03.75. Kk, 03.75. Nt, 05.60. Gg 
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I. INTRODUCTION 



Spinor Bose gases attracted a great deal of attention in the last decade as they exhibit a 
much richer variety of macroscopic quantum phenomena than spinless bosons (see jl| for re- 
view). The physics of three-dimensional spin-1 and spin-2 bosons is rather well investigated, 
both theoretically j2"9| and in experiments with Na and ®^Rb atoms 10-14|. The structure 
of the ground state strongly depends on the interactions, and in particular ferromagnetic, 
polar (singlet-paired), and cyclic phases have been analyzed on the mean field level and 



beyond the mean field The spinor physics of 3D spin-3 bosons is described in Ref. |l5l| 
and, after successful experiments with Bose-Einstein condensates of ^^Cr atoms {S = 3) [l6|, 
experimental studies of the spinor physics in this system are expected in the near future. 

The observation of non-ferromagnetic states requires very low and stable magnetic fields 
(well below 1 mG) at which the interaction energy per particle exceeds the Zeeman energy. 
Presently, the obtained stable field on the level of 0.1 mG is expected to reveal a transition 
between ferromagnetic and non- ferromagnetic states in chromium |l7|, and experiments 
using the magnetic field shielding and aiming at even lower fields are underway 18 1. 

It is important to emphasize that a change of magnetization of an atomic spinor gas under 
variations of the magnetic field requires spin-dipolar collisions, since the short-range atom- 
atom interaction does not change the total spin. In dilute gases of sodium and rubidium the 
spin-dipolar collisions are very weak, and the magnetization does not feel a change in the 
magnetic field on the time scale of the experiment. On the contrary, in a gas of chromium 
atoms which have a large magnetic moment of Qfis, the spin-dipolar collisions efficiently 
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change the magnetization and the gas becomes sensitive to the linear Zeeman effect 

Spinor Bose gases in one dimension (ID) are in many aspects quite different from their 2D 
and 3D counterpats and represent an interesting example of strongly correlated quantum 
fluids. In this paper, having in mind the gas of chromium atoms {S = 3), we assume 
that the system is sensitive to the linear Zeeman effect. We consider a ID spinor gas of 
bosons where the dominant interactions are the density- density and the attractive pairing 
interactions. This choice is justified by the fact that in ID only the latter interaction gives 
rise to a nontrivial quasi-long-range order. In contrast to 2D and 3D, in one dimension 
pairs with nonzero spin S do not condense. This is related to the fact that for S ^ the 
symmetry of the condensate order parameter is non-Abelian. It is well known that strong 



quantum fluctuations in ID dynamically generate spectral gaps for non-Abelian Goldstone 
modes which leads to exponential decay of the correlations (see, for example, j2Q|). As 
far as the polar phase (the condensate of S" = pairs) is concerned, it can be formed 
because the symmetry of the order parameter is Abelian. However, in ID its magnetic 
spectrum is quite different from that in 2D and 3D: in the absence of magnetic field the spin 
excitations have a gap. For a large spin S, the saddle point approximation gives a physically 
transparent description of the polar phase. A sufficiently large magnetic field closes the gap 
and leads to the transition from the singlet-paired (polar) phase to the ferromagnetic state. 
The presence of the spin-gap strongly changes the physics of the ID polar phase and the 
polar-ferromagnetic transition compared to higher dimensions discussed for spin-3 bosons 



in Ref. 



15|. We investigate the ID polar phase and this quantum transition and discuss 



prospects for their observation in chromium experiments. 
II. THE MODEL 

As the atom-atom short-range interaction conserves the total spin, the Hamiltonian of 
binary interactions for (ID) bosons with spin 5* can be written as a sum of projection 
operators on the states with different even spins S of interacting pairs 

2S 

^ = \l 'i^Y.lsM^). (1) 

where x is the coordinate. For the ID regime obtained by tightly confining the motion of 
particles in two directions, the interaction constants 75 are related to the 3D scattering 
lengths a^ niS ) at a given spin S of the colliding pair. Omitting the confinement induced 



resonance 



2l| we have: 



2h: 



2 



where Iq = {h/MuoY^'^ confinement length, M is the atom mass, and uq the confine- 

ment frequency. 

Imagine that all 75 are equal to each other (75 = 7 > 0), except for 7^ at S" = g. 
We then use the relation Yls'^si^) =■ : where fi is the density operator and the 

symbol :: denotes the normal ordering, and reduce the interaction Hamiltonian to the form 
V = (1/2) J dx{'y : n^(x) : +(7^ — ^)Vq{x). For a positive value of (7^ — 7) the system is an 



ordinary Luttinger liquid, but for (7^ — 7) < the situation may change. In 3D a negative 
value of (7g — 7) would lead to a spontaneous symmetry breaking with a formation of the 
order parameter in the form of a condensate of pairs with total spin q. In one dimension 
only a quasi-long-range order is possible and only if g = when the symmetry in question is 
an Abelian one |20|]. Therefore, interactions with negative coupling constants, which have q 
different from zero or from 25" will not produce quasi-long-range-order. The case oi q = 2S 
is exceptional because it corresponds to a ferromagnetic state where the order parameter 
(the total spin) commutes with the Hamiltonian. Therefore, at T = this state can exist 
even in ID. We do not discuss this interesting state, and the only possibility that remains 
is g = 0. So, in our model we have a (repulsive) density- density interaction and the pairing 
interaction that gives rise to the formation of singlet pairs. 

In realistic systems the coupling constants 75 are not equal to each other, although 
they are generally of the same order of magnitude. We thus have to single out the density- 
density interaction in a proper way and then deal with the rest. For example, the interaction 
Hamiltonian ([1]) can be represented as a sum of squares of certain local operators as is usually 
done in the theory of spinor Bose gases [ll. [l5|: 

V = dx Co : h^{x) + ciF^(x) + C2po(a;) + caTrO^x) : +... , (3) 

where F = ip~^Sip, Oij = ip~^{S^S^)ip, the constants Cj are linear combinations of 75, and the 
symbol ... stands for higher-order spin terms which we do not write. The operators : F^(x) : 
and : d^{x) : are given by J^sl^i^ + l)/2 - S{S + l)]Vs{x) and J^si^i^ + l)/2 - S{S + 
l)Y'Ps{x) respectively, where the summation includes all values of S from zero to 25*. We 
then move the S = part of these terms to the term C2Pq{x) and do the same procedure 
with higher order spin terms, which changes the constant C2. The : F^(a;) :, : 0'^{x) : etc. 
terms then no longer contain the interactions with 5 = and, hence, can only lead to 
renormalizations of the density- density and pairing interactions. 

In the case of ^^Cr we have Cq = 0.6576, and the 3D scattering length is ag = 112aB 
where as is the Bohr radius. The exact value of the 3D scattering length a3£)(0) is not 
known and, hence, the constants 70 and C2 are also unknown. In this paper, when discussing 
^^Cr atoms we omit the : F^(x) : and : 0'^{x) : (renormalized) terms, treat C2 as a free 
parameter and focus on the case of C2 < 0. 

We then write down the following Hamiltonian density in terms of the bosonic field 



operators 

j j J j 

where the spin projection j ranges from —S to S, the couphng constant go is assumed to be 
positive, and we put h = 1. The couphng constants go and g are related to cq and C2. For 
example, in the case of ^^Cr we have g = 7co = 4.5576 > and go = — C2. 

III. ZERO MAGNETIC FIELD. SADDLE POINT APPROXIMATION 

We now consider the case ofA^ = 2S' + l >> 1 and apply the 1/A^-approximation to the 
model described by the Hamiltonian density (jl]). First, we decouple the pairing from the 
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density-density interaction by the Hubbard- Stratonovich transformation 

^ [ E ^Z^^] ' ^ ^^V2^ + iA E ^t^,, (5) 
j j 

where A(r, x) and A(r, x) are auxiliary dynamical fields. At large the path integral is dom- 
inated by the field configurations in the vicinity of the saddle point A(r, x) = A, A(r, x) = 
iXo- The values of A and Aq are determined self-consistently from the minimization of the 
free energy. The stability of the saddle point is guaranteed by the fact that the integration 
over the \E', \E'"'' fields yields a term proportional to and therefore the entire action is ~ A^. 
The presence of large in the exponent in the path integral suppresses fiuctuations of the 
fields A and A, thus making the saddle point stable. 
The bosonic action at the saddle point is 

\ (-1)™A+ -iuj-e j \ ^>\^^ ) 

where e = k'^/2M — yU and fi = ^o ~ Aq, with /io being the bare chemical potential. From 
Eq. ([6]) we find the mean field spectrum of quasiparticles (we assume that /i < 0): 



E{k) = ^\/ik' + ^')ik' + kiy, K' = 2M{-fi-A), kl = 2M{-f, + A). (7) 



The saddle point equations are: 



A f dwdfc A 



(9) 



go J {27r)2u^ + E^{k)'- 
diodk e{k) - E{k) 



n 



(27r)2 uj^ + E^{k) ' 
= /^o - gn, (10) 

where n is the density of one of the bosonic species. 

The quasiparticles (spin modes) constitute a (2S'+l)-fold degenerate multiplet. As follows 
from Eq. ([7]), the quasiparticles have a nonzero spectral gap 

E{0) = m = koK/2M. (11) 

This result agrees with the one for = 3 obtained in Ref. |2^]. This is a special feature 
of one dimension. In 2D and 3D the integral in the saddle point equation ([8]) does not 
diverge at small k for /t — )■ 0, and such a gap is not formed. Therefore, one has a gapless 
spectrum of spin modes, which for S* = 3, 2 and 1 has been obtained in the studies of spinor 
Bose gases (see, e.g. and references therein). We would like to emphasize the fact that 
although Eqs. ([8]), ([9]), and (ITOl) resemble the equations for a superconductor, due to the 
bosonic nature of the problem the order parameter amplitude A is not equal to the spectral 
gap, and the latter is related to the parameter k. 

After the integration in Eqs. ([8]) and ([9]) we get the saddle point equations in the para- 
metric form: 



2™ = ^A-ft/l-?i)-*o£(,/l-ij) (12) 



Mgo 




(13) 



where K{x), E{x) are elliptic functions 23|]. From the form of these equations it is clear 



that the ratio k/Zeq is a function of the parameter 
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and kf) can be written in the form: 



(14) 



ko = h{r^){nMgofl\ (15) 




FIG. 1. The ratio k/Zcq as a function of r/. 
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FIG. 2. The parameter & as a function of r/. 



Accordingly, Eq. ( ITTl) for the gap takes the form: 



m = —b {v)-j-{v), 



SO that the gap in units of ngo depends only on the parameter rj. 



m/ngo 
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FIG. 3. The gap m in units of ngQ as a function of r]. 



In the limit of weak interactions where <^ 1, we obtain: 

K 4 /I' 
— ~ - exp — 
ko e \ 7]^ 

and Eq. (ITT]) gives an exponentially small gap: 



m 

e 

For strong interactions, ^ 1, we have 

ko ~ 

and the gap is given by 

m ~ 



exp 



1; b 



V. 

TCT], 



(17) 



(18) 



(19) 



(20) 



The numerically obtained dependence of k/Zcq is displayed in Fig. 1, and the function b{ri) 
is shown in Fig. 2. The gap is presented in Fig. 3. The asymptotic formula f|T8|) obtained 
in the limit of small t] already works with 20% of accuracy for t] = 0.05. With the same 
accuracy the large- 77 asymptotic formula fj20|) is already valid for rj = 1. 

In the limit of weak interactions, taking into account that ^ A and using Eqs. 
and (fT5l) . we get fj, ~ —ngo/AM. Substituting this relation into Eq. (ITOl) we obtain 



n 



9-90 
8 



(21) 



Hence the system is thermodynamically stable for g > Qq. 

The only gapless excitation of the system is the phase mode of the complex scalar field 
A. This excitation describes sound waves of the pair condensate. The effective Hamiltonian 

for the phase mode $ is 

Hphase = \ j dx[KJl^ + K:\d,^f'\, [n(x),$(y)] = -i5(x-y), (22) 

where H is a canonically conjugate momentum. The velocity v and Luttinger parameter Kg 
are extracted from the functional derivatives of the saddle point action and are given by the 
following equations: 

dudk dG{uj, k) dG{-uj, -k) 



2(27r)2 du duj 



iVA2 M3 



/da;d 
2(2^ 



)2 (a;2 + e2-A2)4 Stt J 



■fi{ko/K), (23) 



AuAk dGjco, k) dGj-u, -k) 
2(27r)2 dk dk 



(27r)2^ ' ' (a;2 + e2-A2)4 Stt 7 ^ ' ' 

2NA^ M 

^f2{kolK) (24) 

where G{u,k) is the Green function of the field "i/j, defined as {{'^j{(x!,k)'i/j',{u!,k))) = 
Sjj'G{uj, k), and 



r 8(1 +.-) + (!-. 

^'^^ Jo ^ [{y' + + x')]'/' 

POO 

h{x) = / dyy 
Jo 



2(|/2 + l)2(y2 + ^.2)2 ^ (5/4)(^2 ^ ^2/2 + l/2)2(a;2 - 1)^ 



[(?/2 + l)(y2 + a.2)]7/2 



The functions fi{ko/K) and j^ikoj k) can be expressed in terms of elliptic functions -E(a/1 — n'^/kg) 



and X(Vl - K^/ko), but the expressions are combersome and we do not present them. In 
the limit of weak interactions we have: 

2k ko 
h{ko/n) = f2{ko/K) = — ; ?] < 1 and < ko, (25) 

whereas for strong interactions these functions are given by 

h{ko/K) = —] f2{ko/K,) = g; 77 > 1 and K ~ /cq. (26) 



The velocity v is given by the relation: 



[- 


'f2iko/Ky 








1/2 1 



(27) 



so that for weak interactions using Eqs. ( 125|) . ( |T5l) and 6 ~ 2 we have: 



M 

For strong interactions equations ( |26l) . (|T5l) . and ( fT9l) lead to 



r7<l. (28) 



J^^V; V»l. (29) 

The dependence of v on the parameter t] is displayed in Fig. 4. 
The Luttinger parameter i^,, follows from the relation: 

^s' = Y^[f,-ij ^h{h/K)Uk,/n). (30) 

The scaling dimension of the A field is d = Kg /An and it decreases very rapidly with t], 
which indicates that the mean field approximation works very well. In Fig. 5 we show the 
dependence of d on rj for N = 7. In the regime of weak interactions the scaling dimension 
is exponentially small and it remains significantly smaller than unity even for ~ 2. 

To conclude this part we give a brief summary of the properties of the paired phase. 
With certain modifications, the properties for an arbitrary large spin S are similar to the 
ones for 5 = 1 described in Ref. 2^. Namely, all single particle correlation functions decay 
exponentially. This follows from the fact that the operator always emits a gapped vector 
excitation (Bogolyubov quasiparticle) from the (25' + l)-fold degenerate multiplet (for S = 1 
it is a gapped triplet). Two-particle correlation functions of the operators ^md their 
Hermitean conjugates (no summation assumed) undergo a power law decay. 



IV. MAGNETIC FIELD AND EXACT SOLUTION 

In order to study the influence of the magnetic field on the properties of the singlet- 
paired phase one can also use the saddle point approximation employed in the previous 
section. However, the saddle point equations become too involved. Therefore we resort to 
non-perturbative methods. 
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FIG. 4. The velocity v in the units of y^lngo/M as a function of rj. 
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FIG. 5. The scahng dimension d of the order parameter field A for = 7 as a function of r]. 
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As was demonstrated in Ref. [25|, the model described by the Hamiltonian density (jl]) 
possesses U(l) xO(2S+l) symmetry. Therefore it is reasonable to suggest that the low 
energy sector of this model is described by a combination of the U(l) Gaussian theory and 
the 0(2 S+1) nonlinear sigma (NLa) model. For S=l this was explicitly demonstrated in 



Ref. 



24j . Both the U(l) theory and the sigma model are integrable and the exact solution 



gives access to the low energy sector of the model. At a special ratio of the coupling 



25| that the entire model 



constants one can get even further, since it was demonstrated 
is integrable at a particular ratio of go/g. Below we restrict our consideration to the low 
energy sector where we are not constrained to this particular ratio. 

As we have said, the 0(N) NLcr model is exactly solvable. In the absence of magnetic 
field its excitations are massive particles transforming under the vector representation of the 
0(N) group. This agrees with our result for model (jl]) based on the large N approximation. 
As is always the case for Lorentz invariant integrable models, all the information on the 



thermodynamics is contained in the two-body S-matrix which was found in Ref. [26[. Con- 
sider = 25* + 1 {S integer) and physical particles which have a relativistic-like spectrum 
e{9) = m cosh 9, p{9) = 'y~^msinh6', with mass (gap) m, velocity v, and the total energy 

E = m^coshei. (31) 

i 

For the particles confined in a box of length L with periodic boundary conditions, the Bethe 
Ansatz equations read: 

exp[it;"^mLsinh6'j] = ]^5'o"^(6'i - 9j) Y\ m > 

i^j ai di - + ivr 



V Ail^ - - iTT ^ - \f} -m X'a\' - A^;^ - 2i7r 

\ A« - 9^, + in A« - Ai? + in aI^^ - A^f + 2ivr ' 



n Xas - A^g_/ - ITT _ -p-r Agg - A^^ -ITT 
as-i=l -^"i ~ -^f^z + ITT 63 = 1 Aag - A^s + 

where So is related in the standard way to the integral kernel K{9): 

1 dln^o(^) 
2i^ d^ 

and 



m-K{9), 



K{9) = / e^'^^/"ir(a;)da;/27r, K{uj) 



1 - exp[-2 


u\/iN-2)] 


1 -|- exp[— 


UJ\] 
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The spectral gap m (the particle mass) and velocity v are related to the bare param- 
eters of the model (jlj). For ^ 1 one can use Eq. f|T6l) for the gap, and in the low- 



energy limit the spectrum ([7]) of the spin modes becomes E{k) = \/m? + vH!^ with v = 
{b/2)\JngQ{l + jk'^jM^ so that in the limit of weak interactions we have v ~ ^Jng^^fM. 

Equations (!32|) constitute a system of S + l coupled algebraic equations for the quantities 
6^1, ...Qn, y^x' ■•■Am! \ ■■■AmL ■■■ ■■■Amg^ Integer numbers are eigenvalues of the 5* 
Cartan generators of the group 0(2S+1). From Eq. f l32l) it is obvious that the total energy 
of the system depends on configurations of A's and through them it depends on spin indices 
of constituent particles. Some of the Cartan generators for the problem under consideration 
were constructed in Ref. ( 25|) where it was also shown that the projection of the total spin 



of the system is 

S' = Sn-^mi. (33) 

i 

The Cartan generators commute with the Hamiltonian, and m^'s are integrals of motion. 
Therefore, the magnetic field which couples to through (!33l) . does not violate inte- 
grability. Moreover, once the field is applied, the energies of all eigenstates with mj 7^ 
go up. Thus, at sufficiently low temperature one may consider only eigenstates with no 
A-rapidities, since their energies decrease in the field. When the magnetic field exceeds the 
spectral gap m, the 6'-rapidities start to condense creating a Fermi sea. In the ground state 
the ^-rapidities are distributed over a finite interval (— -B, B\ The distribution function p[&) 
and magnetization are determined by the following integral equations: 

/ K{e - e')pie')de' = — cosh e, (34) 

J-B 27r'f; 

[ K{9 - 9')e{9')d9' = m cosh 9 - hS, e{±B) = (35) 

J~B 

S'/L = S [ d9p{9), (36) 



B 



where h = g^fiB^, with B being the magnetic field and ql the Landee factor. There is 
obviously one transition at = m/S. The magnetization is zero for h < he and it gradually 
increases with the field for h > he (there is another transition in high magnetic fields 
corresponding to the saturation of the magnetization, but the low energy theory cannot 
describe it). In order to find the magnetic field dependence of the magnetization near the 
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transition, where << 1, we approximate the kernel as 

K{e) = 6{e)-A + 0{e^), A = ln4 + ^{l/2 + 2/N -2) -2), (37) 

and look for the solution of Eqs. f l34|) and f l35|l in the form: 

e(^) ^ a{e^ - B^), p{e) ^ const. (38) 
Substituting e{9) and p{9) given by Eq. ( l38l) into Eqs. ( IMl) and ( l35l) ) we get: 



S' _ S 

L 7TV 



2m{hS-m) ^^^|l + 8A/3[2(/i^ - m)/m]^/^ + 0((/i^ - m)/m)}; h > he- (39) 
Keeping only the leading term in Eq. fl39p and restoring the dimensions we have: 



S' V2SmfB-Bc 



L Txhv \ B, 



■'c 



1/2 



B > B,, (40) 



with the critical magnetic field given by 



Equation ( l40l) shows a typical field dependence of the magnetization for the quantum 
commensurate-incommensurate transition. This transition was first studied by Japaridze 



and Nersesyan [27| in the context of spin systems with a gap, where (as in our case) it is 
driven by the magnetic field. Later Pokrovsky and Talapov considered such a transition 
in the charge sector, where it is driven by a change in the chemical potential 28|. The 



magnetization is exactly zero below the critical field and increases as y/B — Be above the 
critical field near the transition. Note that this is quite different from the 3D case where the 
magnetization decreases continuously with the magnetic field when the latter goes below 
the critical value (see, e.g. fisl). 



V. CONCLUSIONS 



We have found that the phase diagram and general properties of ID bosons with a 
large spin S resemble the properties of spin-1 bosons. For the attractive pairing interaction 
[go > 0) the bosons with opposite spin projections create pairs which bose-condense giving 
rise to quasi-long-range order. The saddle point approximation based on the condition of 
large = 25* + 1, gives a transparent picture of the emerging polar phase. The peculiarity 
of one dimension is that all spin excitations have a spectral gap. Hence, in the absence of 

14 



magnetic field there is only one gapless mode corresponding to phase fluctuations of the 
pair quasicondensate. Once the magnetic field exceeds the gap, another quasicondensate 
emerges. This is the condensate of unpaired bosons with spins aligned along the magnetic 
field. The spectrum then acquires two gapless modes corresponding to the singlet-paired 
and spin-aligned unpaired bosons, respectively. At T=0 the corresponding phase transition 
is of the commensurate-incommensurate type, which is qualitatively similar to what we 
have in the case of the 0(3) NLa model. There is a second transition at high magnetic 
fields corresponding to the saturation of the magnetization. However, it is not described by 
the low-energy theory and is beyond the scope of this paper. In the context of ultracold 
quantum gases, the commensurate-incommensurate transition has also been discussed for 
one-dimensional spin-3/2 fermions in the presence of the quadratic Zeeman effect 29|. 



The observation of the commensurate-incommensurate phase transition in a ID gas of 
^^Cr atoms would require (aside from a positive value of go and, hence, a negative C2) fairly 
strong interactions corresponding to the parameter Mgo/n ~ 1, where n = 7n is the total 
density, so that rj ~ 0.5. Then the spin gap in the polar phase is of the order of ng^ and 
(assuming gQ by a factor of 3 smaller than g) can be made on the level of 100 nK at ID 
densities n ~ 10^ cm~^. Then the transition occurs at the critical field of the order of 0.2 mG 
and can be observed at temperatures T ~ 20 nK. This however is likely to require the ID 
regime with a rather strong confinement in the transverse directions (with a frequency of the 
order of 100 kHz as in the ongoing chromium experiment in the ID regime at Villetaneuse 



17|). 



Note added 



After this work has been finished the Villetaneuse group has reported the observation 
of the demagnetization transition for ^^Cr atoms in the ID regime, under a decrease of the 
magnetic field to below 0.5 mG. However, the experiment is done at temperatures ~ 100 nK 
and the state which is reached by decreasing the magnetic field does not necessirily reveal 
the nature of the ground state due to thermal excitations (and due to diabaticity at the 
transition in the experiment). 
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